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LEWICE Ice Accretion Prediction 


LEWICE is a software package the predicts the size, shape, and location 
of ice growth on aircraft surfaces exposed to a wide range of icing 
conditions. 


LEWICE also models the 


¥ ___ Flow solution using potential flow or structured viscous behavior of thermal ice 
solver protection systems while 
v Particle trajectory calculation, including impingement exposed to the same range 
limit search for collection efficiency and multiple drop — a 
size distributions of icing conditions. 
v Integral boundary layer routine calculates heat transfer 
coefficient 
S@ar 
v Quasi-steady analysis of control volume mass and “ Re cA 


energy balance in time stepping routine 

v Geometry modification using density correlations to 
convert ice growth mass into volume allows multiple 
time-step solutions 

v All physical effects modeled, including turbulence, 
bouyancy, droplet deformation, breakup and splashing 

v Extensive validation against experimental data 
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Adhesion Models 
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Penn State Tests 


Most shedding models my roy 
used today are based on J 
empirical adhesion data, 
and are mainly a function 
of temperature... 


SAE 2007-01-3315 


NASA TM-107162 


TOTAL TEMPERATURE. 


SAE 2003-01- 2093 


...but need to account for geometric conditions as well. 
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Typical Shedding Model 


= _ 2 
Fe _ F apecion - Fe ieion _ Mj-.¥G2 


entrifugal 


F Adhesion = Or iescs 


Fe Cohesion a og en eee 


' Cohesive tensile stress 


Ice fragment 


Fortin, G., and Perron, J.; “Spinning Rotor Blade Tests in Icing 
Wind Tunnel.” AIAA 2009-4260, 1st AIAA Atmospheric and Space 
Environments Conference, June 22-25, 2009, San Antonio, TX. 


Itagaki, K.; “Mechanical Ice Release Processes — I. Self-Shedding from 
High Speed Rotors.” CRREL Report 83-26, U.S. Army Cold Regions 
Research and Engineering Laboratory, Hanover, NH, October 1983. 
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NASA Glenn Icing Research Tunnel 
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Rotor Model in the Icing Research Tunnel 
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Rotor Model in the Icing Research Tunnel 


Direction of Air Flow > 


— = 
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Workflow of the STAT Script 


Define Parameters 


¢Frames per second/timestep 
¢Blade rotation per timestep 
e4 points for image transform 
eRotor disk diameter 


Subtract Background 


eUsing median image of all shed 
images 


| Indentify Shed Images 


| elmages masked to remove rotor disk 


Peaks detected in image intensity 
elmages Imported 


Estimate Shed Origination 


eFind tangent lines to shed particles 


eldentify shed angle at first tangent 
point 


Images Processed 


¢Masking (remove rotor disk) 

eMedian filtered 

eOtsu threshold to create binary 
image 

eMask within 15 pixels of Otsu 
thresholded gradient image 


Transform Images 


eUsing 4-point projective plane 
transformation 


Identify Particle Statistics 


eArea determined by pixel count 
eCentroid located 
eVolume estimated from area 


Calculate STM Fronts 


¢User input to fit model 
eSpecify break points on leading edge 


Plot Output 


ePlot binarized image 
ePlot STM shed fronts 
¢Plot rotor disk line 


¢Plot background subtracted rotor 
disk 
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Finding Shed Events 


Description Rotor Speed Blade Pitch Tunnel Static Temp. Tunnel Speed LWC MVD 
(rpm) CEs) (°F) Gs) (gim*)__(um) 

Chordwise De-ice 1200 10 60 0.5 15 
Spanwise De-ice 1200 5 60 0.5 15 
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Image Processing — Original (Run 71) 


X7 


Y watnal = 60 kts 
RPM = 1200 
Collective = 5° 
Tytatic = -4 °F 


*static 
= 6? De 
spray wa 2) 
LWC = 0.5 g/m? 
MVD = 15 micron 


Frame Rate = 48 


Spanwise De-Ice Evaluation Run 
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Image Processing - Transformed 
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Image Processing — Background Subtracted 
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Image Processing — Background Subtracted 
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Image Processing — Final Stage 
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Run 67 Shed Event 


Chordwise De-Ice Evaluation Run 
V — €/ kte 

V tunnel — 60 kts 
RPM = 1200 
Collective = 10° 
T =. °F 
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Frame Rate = 320 {ps 
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Run 67 Shed Frame -4379 
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SHEDDING TRAJECTORY 
MODEL (STM) 
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Terms (1/2) 


L =rotor disc radius 
d, = start of shed ice 
d, = end of shed ice 
X, = /adial position of shed ice, inner edge 


8 =rotational position of the rotor, 
— Typically 6; = 0 at shed start 


— §,can be adjusted to correct STAT’s automatically 
identified shed start 


Z =L-—d,+d, 


— When ice reaches blade tip and breaks off 


Aw 
Shown on next slide: J 


A =cross-sectional area of ice, 
normal to X 


W =width of A at interface 
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Terms (2/2) 


e Ais the blue area 
e¢ Wis the length of the red line 


e Mis the mass of the ice 
e F, is the centrifugal force on the ice 


e F,, are the non-conservative forces 
— bis coefficient for viscous damping 
— F; is coefficient for friction (per unit length) 


: p= et re 
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On-Rotor Position Calculation 


Centroid Interfacial 
Mass Position Area 


So 


ia 
pA(d> — dy) ( + a) w* — (dz — d,)(Wbx, + Fy"), x5 <z 
DiAxs) = 


Ll? -—x 
a we —-(L- cy i “ tak ite = 
baa re a oer 
: _ Ti Ti Ti 
M; (xs) we _ ES — erm erm erm 
w* (x,bW + F;) 
(20 (ds = 4) ) = 7 
a;(x;) = Coe ee) pA : 
T°" Mi(xs) w? (x,bW + Fr) 
7 ye) any Cds = 
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On-Rotor Position Calculation 


The on-rotor position is necessary to get 
initial velocity of ice leaving the rotor 


Problem: 
2 
Ww (d, — ad, = 2F) 
X. + BX. — @*x. = —_————,x%. <2 
: : 2 : (—B + VB? + 4@2) 
, a w*(L — 2F) "12 = 2 
X, + BX, -—x, = ———_,%, = Zz 
2 2 (—B + VB? + 2w?) 
a ar 
: _ d, +d, —2F 
Solution: ~~ 9a fy 
2 
rer eae 3 
ds =o i Spat eee 
pe Ge ee Sr 77 ee : 21-2 
.(t) = 
c3e73§ + c,e™! + 2F —L,x, =z a ull ak aa am Tae) —— — 4 
3°” (I 4 
Vig mirgl2l —2F = ay + dj) 
—— iTS 
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Off-Rotor Position Calculation 


x; = Lcos(wt + 6;)i+ Lsin(wt + 0;) 7 
Vt = wL 


Cine ones <2 


== 
: a Ne ae 


C3P3e 


V;(B, F) = V(B, F) cos(wt + 6;) — V; sin(wt + 6;))i 
+(V,(B, F) sin(wt + 6;) + V; cos(wt + 6;))f 


Predicted position of shed ice with quadratic drag 


(c is the drag term): 
2m 


m cV;(B, F)t 
xX>(t,B, F,c) = —In|{ 1 + ————— ] + x; 
2 P a i 
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RESULTS 
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Estimated Shed Ice Mass for High Speed Impact 
Analysis 


Initial ice shape Projected area to diameter 
Assuming , ndéq 
0 4 A= Nea — 
0 OO f 
—+> doooo000 ——> 4, > 3 
goooo0o0 4 (Ay 
oo0000 ea) 


jE 
shed (kg) | particles | mass (kg) | masses (kg) 
manell inc | sh.) atits. | saa 
(chordwise) 0.0511 38 0.00134 0.0023 


fepame)| orese | sos | oem | oy 

(spanwise) 0.1064 108 0.00099 0.0037 

jiraory | ooos | «| ooous | ons 
(natural) 0.0045 4 0.00113 0.0006 
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Estimated Mass Distribution 
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Conclusions 


¢ Two runs consisting of three shed events were 
analyzed 


¢ The STAT tool was developed to post-process the 
data 


¢ A model (STM) to predict the trajectory of a shed 
front was developed 


— Viscous forces at the interface between the ice and the rotor 
are significant 


— Bis the most important fit parameter, B = 60 provides 
current best fit 


— There is insufficient data to estimate the importance of the c 
term 


— Interstitial breaks increase the shed fan angle 
¢ The mass distribution of shed ice was estimated 
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Questions? 
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